Convergence of the Linear 5 Expansion 
in the Critical O(N) Field Theory 
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The linear 8 expansion is applied to the 3-dimensional O(N) scalar field theory at its critical point 
in a way that is compatible with the large- N limit. For a range of the arbitrary mass parameter, 
the linear 8 expansion for (4> 2 ) converges, with errors decreasing like a power of the order n in 8. 
If the principal of minimal sensitivity is used to optimize the convergence rate, the errors seem to 
decrease exponentially with n. 
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Very few systematically improvable methods are 
available for calculating nonperturbative quantities in 
field theory. Such a method has the advantages that 
more accurate results can be obtained at the expense of 
additional effort and that reliable error estimates can be 
made. One systematically improvable nonperturbative 
method that is well-developed is the Monte Carlo method 
applied to the field theory formulated on a discrete lat- 
tice. Since this method is not universally applicable, it 
is important to develop other nonperturbative methods. 

A classic nonperturbative problem that was recently 
solved definitively is the shift due to interactions in the 
critical temperature T c for Bose-Einstein condensation 
(BEC). If the potential between two bosons is short- 
ranged, the leading order shift is linear in the s-wave 
scattering length a: AT C /T C = cn x / z a, where n is the 
number density of the bosons and c is a numerical con- 
stant. Baym et al. Q showed that the coefficient c could 
be determined by a nonperturbative calculation at the 
critical point of an effective 3-dimcnsional statistical field 
theory with 0(2) symmetry. Lattice Monte Carlo calcu- 
lations by Kashurnikov, Prokof'ev, and Svistunov and 
by Arnold and Moore [|| give the result c = 1.32 ± 0.02. 
The second order correction to AT C /T C proportional to 
(an 1 / 3 ) 2 has also been calculated J|]. The definitive solu- 
tion to this problem makes it useful as a testing ground 
for other nonperturbative methods. 

One systematically improvable nonperturbative 
method that has been applied to the shift in T c is the 1 /N 
expansion. The 0(2) field theory relevant to BEC can be 
generalized to O(N). Baym, Blaizot and Zinn- Justin cal- 
culated the coefficient c analytically in the large- N limit 
and obtained c = 2.33 The first correction in the 
1/N expansion reduces c to 1.72 ||. These results seem 
to be converging to the lattice Monte Carlo result. The 
analytic result in the large-TV limit can be used to test 
other nonperturbative methods. If a method fails to give 
the correct answer in the large- N limit, one should be 
suspicious of its predictions for N = 2. 

Another nonperturbative method that has been ap- 



plied to this problem is the linear 5 expansion ||, also 
known as optimized perturbation theory J7J or variational 
perturbation theory . In this method, an arbitrary pa- 
rameter m is introduced into the theory and calculations 
are carried out using perturbation theory in a formal ex- 
pansion parameter 8 = 1. The rate of convergence can 
be improved by adjusting m at each order in 6 using the 
principal of minimal sensitivity (PMS). It was first ap- 
plied to the calculation of AT C by de Souza Cruz, Pinto 
and Ramos §. At 2 nd , 3 rd , and 4 th order in S, they ob- 
tained c = 3.06, 2.45, and 1.48, respectively H, which 
seem to be converging to the lattice Monte Carlo result. 

The fact that corrections can be calculated system- 
atically is no guarantee that they actually improve the 
result. Hopes for the convergence of the LDE are based 
largely on studies of its application to the anharmonic os- 
cillator. The LDE has been proven to converge for order- 
dependent choices of m that include the PMS criterion as 
a special case. The finite temperature partition function 
converges exponentially, with the errors at n th order de- 
creasing as exp(— bTn 2 / 3 / g 1 / 3 ) where g is the strength of 
the anharmonic term in the potential and b is a numeri- 
cal constant |l(J. The energy eigenvalues E n have been 
proven to converge uniformly in g as n — ► oo fllf . In par- 
ticular, the leading term in the strong-coupling expansion 
for the ground state energy Eq converges exponentially, 
with the errors decreasing like exp(— b'n 1 ^ 3 ) |p^| . 

The anharmonic oscillator is equivalent to a Eu- 
clidean field theory with a single real-valued field in 1 
space dimension. The statistical field theory relevant to 
BEC is a generalization to a multicomponent field in 3 
space dimensions. The most serious obstacles to general- 
izing the convergence proofs for the anharmonic oscillator 
to this more complicated problem come from the infrared 
(IR) and ultraviolet (UV) regions of momentum space. 
It is reasonable to expect the convergence behavior to be 
similar if appropriate IR and UV cutoffs are imposed on 
the field theory. The coefficient c is insensitive to the 
UV region, so we do not expect any complications from 
taking the UV cutoff to oo. However, c is very sensitive 



to the IR region, so convergence in the presence of an 
IR cutoff gives no information about the behavior of the 
LDE in the limit as the IR cutoff goes to 0. 

In this letter, we present evidence that the LDE is 
indeed a systematically improvable method. We show 
that it can be implemented in a way that is compatible 
with the large- TV limit. We present evidence that the co- 
efficient c in the large- TV limit converges to the analytic 
result of Ref. for a range of to, with errors that de- 
crease as a power of the order n in the LDE. If the PMS 
criterion is used to optimize the convergence rate, the 
errors seem to decrease exponentially in n. 

The lagrangian density for the O(N) field theory 
relevant to BEC is 
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where is the self-energy and f p = J d 3 p/(2,Tr) 3 . The 
integral over p is divergent and requires a UV cutoff. The 
critical point can be reached by tuning the parameter r 
to the value r = — S(0). This condition reduces to r = 
if u = 0. The difference A between the critical values of 
(d 2 ) at a nonzero value of u and at u = is 



A = TV 
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The integral over p is convergent and therefore no longer 
requires a UV cutoff. At the critical point, the only rel- 
evant length scale is set by the parameter u. Since A 
has dimensions of length, it must be proportional to u 
by dimensional analysis. The determination of the coef- 
ficient of u requires a nonperturbative calculation. In the 
large- TV limit defined by TV — > oo, u — > with TV it fixed, 
the coefficient is known analytically [Q: 



A = -TVu/(96tt 2 ) (large TV) . 



(4) 



The coefficient in the expression for the shift in T c is 
c = -I28tt 3 C(§)~ 4/3 A/u, with A evaluated at TV = 2. 

The linear 5 expansion (LDE) is generated by a la- 
grangian whose coefficients are linear in a formal expan- 
sion parameter 5: Cg = (1 — 6)£q + 6C, where Co is 
the lagrangian for an exactly solvable theory. The la- 
grangian Cs interpolates between Cq when <5 = and C 
when 5 = 1. To apply the LDE to the O(N) statistical 
field theory defined by the lagrangian (Q), we choose the 
solvable field theory to be the free field theory with mass 
to. The lagrangian can be written Cs = Cq + C- m t, where 
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FIG. 1: The 4 th in the series of diagrams for E(p) that survive 
in the large- TV limit. 



Calculations are carried out by expanding in powers of 
5, truncating at n th order, and then setting 6=1. 

In order to apply the LDE to the shift in T c) we need 
a prescription for generalizing the quantity A defined in 
(H) to the field theory defined by the lagrangian Cs ■ The 
prescription must have a well-defined expansion in pow- 
ers of 5, and it must reduce to (||) when 5=1. The 
simplest prescription is to use the expression (|^), where 

is the self-energy for the field theory with lagrangian 
Cs- In the previous application of the LDE to the shift in 
T c , the authors used an alternative prescription with the 
additional term to 2 (I — 5) in the denominator of the first 
term in the integrand of (||) || . This prescription has the 
correct limit as 5 — * 1. It can be expressed as an integral 
with a well-defined large- TV limit plus an additional term 
— (N/4:ir)my/l — 5. Because of this additional term, the 
limit 5 — ► I does not commute with the large-TV limit. 
Using this prescription, the prediction for A defined by 
the PMS criterion at n th order in the LDE, scales like 
jy2—\/n u ag _> qq ^ w hile the correct result is TVu. 

The prescription (|^) defines A (it, m, 5) as a function 
of three variables. The n th -order approximation in the 
LDE is obtained by truncating the expansion in 5 at n th 
order to obtain a function A'™) (u, to, 5) and then setting 
5 = 1. At any finite order n, the prediction of the LDE 
depends on to. As to varies over its physical range from 
to +oo, the range of the prediction for A extends out 
to ±oo depending on the order in 5. Some prescription 
for to is required to obtain a definite prediction. The 
PMS prescription is (d/dm)A^ (u, to, 5 = 1) = 0. After 
setting 5=1, A^"^ is a function of u and to only. By 
dimensional analysis, the value of A*-™) at a solution to to 
the PMS criterion is proportional to u. By allowing the 
variable to to change with the order n, the PMS criterion 
may improve the convergence rate of the LDE. 

Using the prescription (||), the LDE for A in the 
large-TV limit can be calculated to all orders in 5. The 
leading contribution to X(p) — S(0) comes from the se- 
ries of diagrams whose 4 th member is shown in Fig. [I]. 
Since S(p) — S(0) is of order 1/TV, the leading contri- 
bution at large N is obtained by expanding (||) to first 
order in E(p) — 11(0). The expression for the large- TV 
diagram for A with n + 1 loops can be reduced to a 1- 



dimensional integral multiplied by to ( Tl l ^u n . In addi- 
tion to the diagrams for E(p) generated by the interaction 
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FIG. 2: A/(-7Vfi/967r 2 ) in the large-iV limit as a function of 
fj, at n th order in the LDE for n = 3, 4, 5, 7, 11, 19, 35. The 
curves for n = 7, 11, and 19 appear in order between those 
labelled n = 5 and 35. 

term 6 u(tfi 2 ) 2 in (|^), we must also take into account in- 
sertions of Sr and — Sm 2 . The effect of the Sr insertions 
is to replace each by — S(0). The effect of the 
—Sm 2 insertions is to replace m 2 by m 2 (l — <5). Summing 
all the large- TV diagrams, we obtain 

/ \ 71— X 

A = 8Nu/ {24n 3 )J2 (-<*/( vT^M)) 

n=2 

/>oo 

x / dyy 2 [A{y)] n - 1 /{y 2 + l) 2 , (7) 
Jo 

where /u = 48irm/(Nu) and A(y) = (2/y) axctan(y/2). 
The prediction for A at n th order in the LDE is obtained 
by expanding (Q) as a power series in <5, truncating after 
order 5 n , and then setting 8 = 1. 

It is easy to show that if the LDE for (pt) converges, 
it converges to the correct analytic result (ffl). After in- 
terchanging the order of the sum and the integral, the 
sum can be evaluated. Upon taking the limit 5 — ► 1, 
all dependence on [i disappears. Evaluating the integral 
over y gives the result (|J). Thus, if the LDE converges 
at some value of fi, it should converge to (||). 

The manipulations that showed the convergence to 
involved several interchanges of limits. It is difficult to 
translate the conditions for the validity of each step into 
a condition for the convergence of the LDE. However, the 
convergence can be easily studied numerically. In Fig. 
we show A as a function of \i for several orders in the 
LDE: n = 3 and n = 2 j + 3, j = 0, 5. The horizontal 
line is the analytic result (|4|). The results are consistent 
with convergence to (|4|) for all \i greater than a critical 
value fi c which we estimate to be /i c rs 0.71. If /i < fi c , A 
seems to diverge to +oo for n even and to — oo for n odd. 
For any fixed /i > /z c , the convergence with n is very slow. 
In Fig. ||, we show a log-log plot of the fractional error 
e n as a function of n. The squares lie close to a straight 



FIG. 3: Log-log plot of the fractional error e n as a function of 
the order n in the LDE. The squares, diamonds, and circles 
are for /i = 1.039, the real solution to the PMS criterion, and 
the solution with maximal |ImA|, respectively. The lines are 
simple curves that pass through the last 2 points. 



line, indicating that the errors decrease like a power of n. 
The dotted line that goes through the last two points is 
e n = 0.70n -0 - 37 . The errors decrease roughly like n -1 / 3 . 

The rate of convergence can be improved by using 
the PMS criterion to choose a value of [i that depends 
on the order in the LDE. At n th order, this criterion is 
a polynomial equation in /i of order n — 2. For n even, 
there are no real roots. For n odd, there is always one 
real root that corresponds to the maxima of the curves 
in Fig. ^. The resulting fractional errors are shown as 
a function of n in Fig. ^. The diamonds lie close to 
a straight line, indicating that the errors decrease like 
a power of n. The dashed line that goes through the 
last two points is s n — 0.78n _0 46 . The errors decrease 
roughly as nT x l 2 . 

Although the PMS criterion at n th order has at most 
one real solution, there are always n — 2 complex- valued 
solutions. Studies of the anharmonic oscillator have re- 
vealed that there are families of complex solutions with 
much better convergence properties than families of real 
solutions jn|. In our problem, at any odd order n, the 
real solution always gives the value of ReA that is far- 
thest from the correct result Thus this family of so- 
lutions gives the slowest possible convergence rate. How- 
ever there is a strong anticorrelation between the errors 
in ReA and ImA. This is illustrated in Fig. ||, which 
is a scatter plot of |ImA| vs. ReA for the solutions to 
the PMS criterion for n — 35. The solutions with the 
most accurate values for ReA are those with the largest 
values for |ImA|. Thus we can define a nearly optimal 
family of solutions by choosing those with the maximal 
values of |ImA|. While ImA for these solutions shows no 
sign of converging to 0, there is dramatic improvement in 
the convergence of ReA. A log-log plot of the fractional 
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Re A I(-Nu/96k ) 

FIG. 4: Scatter plot of |ImA| vs. ReA for the solutions n of 
the PMS criterion at 35 th order in the LDE. 



errors in ReA is shown in Fig. |^. The downward curva- 
ture of the points indicates that the errors decrease faster 
than any power of n. The solid line that goes through 
the last two points is e n = 0.32(0.93)". The errors are 
decreasing faster than this exponential. 

We have calculated A for finite N to 4 th order in the 

Setting N 



LDE using our prescription 1 14 . Setting N — 2 and using 
the PMS criterion, we obtain a real result A^ 3 ) = 0.192 
at 3 rd order and a complex result = 0.214±0.084i at 
4 th order. Taking A' 3 ) and ReA^ 4 ^ as the predictions, the 
corresponding values for the coefficient in AT C are c = 
0.447 and 0.492. They seem to be slowly approaching the 
lattice Monte Carlo result c = 1.32 ±0.02 [f| from below. 
The errors in the 3 rd and 4 th order results are 66% and 
62%, which is larger than the errors of 52% and 44% in 
the corresponding predictions in the large- N limit. The 
percentage improvement in going from 3 rd to 4 th order 
is only half as large as in the large- N limit. The 3 rd and 
4 th order predictions for N = 2 using the prescription 
of Ref. H have errors of 85% and 20%. The small error 
in their 4 th order prediction may be fortuitous. The low 
order predictions using their prescription are skewed by 
the term in A that does not have a well-behaved large- iV 
limit. The prescription for A obtained by deleting that 
term is equally valid, and it gives values below the lattice 
Monte Carlo result with errors of 68% and 65%. 

Although the convergence rate of the optimized LDE 
appears to be exponential in the large- N limit, it is still 
rather slow. One must calculate to about 18 th order in 
5 to achieve 10% accuracy. For general N, it may be 
feasible to calculate A to 5 th order, but it would be diffi- 
cult to go to much higher order. The slow convergence in 



the large- TV limit suggests that even if the LDE also con- 
verges for N = 2, a strict expansion in 5 is not useful for 
quantitative calculations. It may however be possible to 
use order-dependent mappings to change the expansion 
in S into a more rapidly converging expansion p"5| . 

In conclusion, we have shown that the LDE for the 
quantity A in the 3-dimensional critical O(N) field the- 
ory in the large N limit converges if we use an appropriate 
prescription. If the PMS criterion is used to optimize the 
convergence, the errors seem to decrease exponentially in 
the order of the LDE. This provides some hope that the 
LDE may also be a systematically improvable approxi- 
mation scheme in the case N = 2 relevant to BEC. With 
the use of order-dependent mappings to accelerate the 
convergence, it may be possible to develop the LDE into 
a general and powerful tool for quantitative calculations 
in superrenormalizable field theories, even at a critical 
point. 
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